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Abstract 

Testing symptomatic individuals for a disease can deliver treatment resources, if tests’ results 
turn positive, which speeds up their treatment and might also decrease individuals’ contacts 
to other ones. An imperfect test, however, might incorrectly consider susceptible individuals 
to be infected (false positives). In this case, testing reduces the epidemic in the expense of 
potentially misclassifying individuals. We present a mathematical model that describes the 
dynamics of an infectious disease and its testing. Susceptible individuals turn to “susceptible 
but deemed infected” at rate 6. Infected individuals go to a state “infected and tested 
positive” at rate a. Both rates are functions of test’s sensitivity and specificity. Analysis of 
the model permits us to derive an expression for Rq and to find the conditions for reaching 
Rq < 1, i.e., when the disease-free equilibrium is stable. We also present numerical results 
to cover interesting scenarios such as using different tests and to compare these results. 
We find for different sensitivity and specificity values the conditions permitting to get the 
basic reproduction number Rq < 1, when originally, i.e., without testing, we would have 
Ro > 1. We also find for a given sensitivity and specificity, the critical testing rate for 
reaching Rq < 1. 


1 Introduction 


Treatment for communicable diseases can be quite costly. In the case of multiresistant tuber¬ 
culosis infection the treatment per patient hast a cost that is much higher than an infection 
from non-resistant strain of the pathogen (Cohen et al. 2009). This kind of treatment, 


however, relies on testing diagnostics which have limitations. They basically might be more 
accurate at the expense of time taken to have the results. It is often best to start treatment 
as soon as the case is diagnosed, but the test results might suffer from the test imperfec¬ 
tions. A model that indicates the cost of an epidemic as a function of the test sensitivity 
and specificity can help mitigate or even avoid an epidemic. 
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The test imperfections appear from the probabilities of finding true positives or true 
negatives. The first effect comes from the specificity, i.e., the probability that the test 
indicates a negative diagnosis given that the tested individual does not have the pathogen. 
If the specificity is not so high a significant fraction might be considered infected, whereas 
those individuals are not. In some cases, starting treatment is the priority over waiting for 
more accurate results. Therefore imperfections here might increase the cost. 

Also, for real infected individuals, it is important that the test has good sensitivity, i.e., 
its rate of positive diagnosis in the case of positive tested individuals should be high. Again, 
this comes at an expense of rapid diagnosis. The effect is that individuals that are correctly 
diagnosed start treatment and reduce both the time they could infect other people and the 
contacts to susceptible individuals, because of isolation and other measures. 

In the case of tuberculosis (TB), for instance, an example is the Xpert (R) MTB/Rif 
assay test for diagnosis of TB developed for the GeneXpert platform (WHO 2011). It 


consists of a molecular testing that can be done on demand and closer to individuals in 
need. A meta-analysis of these reports was done by Walusimbi et al. (2013). The test 
results take generally about three hours, whereas a more established protocol (sequence of 


tests and surveys) would result in treatment starting only after a few days. Walusimbi et al. 
(2013) found pooled specificity and sensitivity estimated at 0.67 (Cl: 0.62 - 0.71) and 0.98 
(Cl: 0.97-0.99), respectively. 

In the case of hepatitis C, de Paula Scalioni et al. (2014) report results for rapid tests 
which can have high values for both specificity and sensitivity, but can vary depending on 
the kind of sample (serum, whole blood, saliva). 


Eaton and Hallett (2014) estimate the effect of early-stage transmission of HIV, he., 


before testing individuals for HIV infection, a stage at which they did not started receiving 
treatment. 

The recent literature on disease modeling and diagnosis is extensive. We point to the 


review by 

Zwerling et al. (2015 

on modeling of TB, investigating diagnosis and therapies. 

In particular, Raimundo et al. 

2014) have a an approach on modeling focused on resistant 

TB. 

3alje et al. 

(2014) discuss t 

re implementation of testing for multiresistant TB in India. 

Also, 

Cox 

(2014 

) in a recent editorial described the need for simple models for the sake of 


interpretation of results, especially for people from outside of modeling. In our model, we 
consider the minimum number of variables to give important insights. 

In this paper we present a modeling approach that permits us to obtain Rq as a function 
of testing specificity and testing sensitivity. Our model is described as a set of ordinary 
differential equations, close to a SIR model. We present analytical results from the model 
and also some numerical simulations that show the cost increase due to testing along with 
the desired effect of the mitigiation or avoidance of epidemics. 


2 Model 

Our ODE model, leveraged from a classical SIR model, has two components: S m , to describe 
the number of susceptible individuals that take the test and are deemed incorrectly to be 
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infected, and I m , that describes the number of infected individuals that are test-diagnosed 
positive and start treatment. As usual, we describe the other compartments: S, to describe 
number of suscetibles, I for number of infected individuals and R for recovered ones. 

Individuals that are susceptible get infected at rate /3. These individuals are called the 
regular susceptibles. The model describes the situation in which susceptible individuals take 
a test and are deemed incorrectly to be infected at rate f3 m . This class is for the “susceptible 
but deemed infected”. We generally assume that /3 > f3 m , for infected ones might be under 
treatment, possibly in isolation, therefore these individuals are less likely to transmit the 
pathogen to another individual. 

Individuals leave the population at rate ju, equal for all classes, because we assume a 
disease that does not cause significantly higher mortality rate. The rate of new individuals 
entering the population ( e.g ., typically by birth) is also given by /r to result in a constant 
population size. We normalize all variables by the population size. 

Individuals pass from susceptible to a susceptible-deeemed-infected state at rate 6. This 
will depend on the rate r at which individuals are tested and the test specificity e. The test 
specificity is given by the probability of true negatives, i.e. P(test- | uninfected) (Newman 
2003). Hence, 8 is the product of the rate r and P(test+ j uninfected) — ! — e, such that 


8 = r(l — e). 

The test sensitivity ip is given by the rate of true positives, i.e. P(test+ | infected) 
(Newman 2003). Also, infected individuals can be tested and, if the test is positive, start 
treatment. Therefore, infected individuals will pass from class / of regular infected indi¬ 
viduals to the class of I m infected under treatment at rate a = rip. Finally, we assume 
that individuals that are under treatment get to recover faster. Therefore, regular infected 
individuals recover at a rate 7 , whereas infected-under-treatment individuals recover at rate 
7 m- Typically, > 7 . 

We present the ODE system that describes this model: 

7 c* 

= -OS - PS(I + I m ) - vS + n 


dt 

dSm 

dt 

dl 


= 6S-p m S m (I + I m )-iiS n 


—j~ — PS(I + Im ) + PmSm{I + Im) — (ct + 7 + AO-f 


(1) 


dirt 


— Oil 'Jmlm f^In 


dt 
dR 

7 , — 7 mlm T 7 1 fJtR 
dt 

Table [l] contains a summary of the variables and parameters, and their respective de¬ 
scriptions as well as values used in this work for investigation of numerical scenarios. 


Cost of infection 

We apply here a model in which treatment cost varies linearly with the number I of infected 
persons and the number I m , but also the number S m of susceptible but deemed infected. 
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Table 1. Model parameters, variables and respective descriptions. 


Parameters 

Description 

Value 

P 

transmission rate in state S 

.15 

7 

recovery rate in state 7 

.1 

7 m 

recovery rate in state I m 

.15 

Pm 

transmission rate in state S m 

.1 

V 

mortality/birth rate 

0.003 

r 

testing rate 

different values 


test sensitivity 

different values 

e 

test specificity 

different values 

a = rip 

rate at which individuals enter state I m 

different values 

6 = r( 1 — e) 

rate for entering state S m 

different values 

Variable 



S 

number of susceptible individuals 

normalized 

Sm 

number of susceptible individuals but deemed infected 

normalized 

i 

number of infected individuals 

normalized 

Im 

number of tested-positive infected individuals 

normalized 

R 

number of recovered individuals 

normalized 


The cost is a function of the rate r and time r, given by C a (r) = JJ wsS m (r, t ) + wjl ( r , t ) + 
w M Im(r,t)dt. The relative cost C(t) is given by the ratio between the cost under rate r 
and the cost at no treatment, i.e., r = 0, as follows: 

, . _ fo w s S m (r , t) + wil(r, t ) + w M Im(r , t)dt 
Jo wsS m (0,t) + 10/7(0,t) + w M Im(0,t)dt ’ 

where res, uq, and tom are weights to each of the variables I, and I m , respectively. 

3 Results 

System dynamics 

We start the analysis by considering that there are only susceptible individuals, without any 
infectious individuals. In this case, we have clearly 7(f) = 0, I m (t) = 0 and the left hand 
side of equations in [l] become zero. The solution is given by 

(Sf , S mJ ,If, I m ,fi Rf) = (— _|_ q ■, ^ q 10,0, 0 ). 

As noted above, he system has an equilibrium in which S m = jJJe- Therefore, a fraction 
of the population demands treatment resources even when there are no infected individuals. 
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We follow from this initial result to derive a matrix that describes the “flows” into in¬ 
fectious states and out from infectious states. Note that the state S m is not really an 
infectious state, since the individual is still susceptible. We apply the method of analyzing 
the Next-Generation Matrix (NGM, described by Diekmann et al. (2009)) K and find: 


K = 


Pr | PmQ 

/ i .-|-0 

Q+7+/Z 


+ 


a 


Pn 


PmO ’N 

M+0 / 


(“+7+M)(7 m +m) 

0 


0^ | PmQ 
7m 

0 


The dominant eigenvalue of the matrix K given by the spectral radius of matrix K 
results in the basic reproduction number Rq, the number of individuals that get infected 
upon a single infected individual: 


Ro 


P^ + Pm@ 


(1 + 


a 


(0 + n)(l + H + ay 7 m + V 


( 3 ) 


Endemic state 


Analysis of the system dynamics permits us to show that it reaches equilibrium in en endemic 
state, as expected depending on the threshold given by Rq. The endemic state is given by: 


S* = 


Q* _ 

m 


I* = 


R* = 




0 + h + + Ck/( 7m + A*)) 


0 


Pm,I* (1 + Oi/{ 7 m + fl)) + flj \0 + fl + PI*{1 + a/('y m + /i)) 
cd* 


7m + At 

——-1- 'y 

7m+/£ ' ^ 

h 


7m 


Using the equations above we can find an equation containing only the number /* of 
infected in the endemic state. Solving the equation for I* requires finding the solutions to 
the polynomial 


Pii 1 ) = 


(0 + [i)fi 


PPm(l + Ct/(7m + tp) 2 


(R 0 - 1) + I 




(0 + /l + fin/Pm) 
'y + a + fi P(1 + a/pjm + //)) 


-J 2 (4) 


Clearly, i?o > 1 is a sufficient condition for a biologically feasible solution (a nonnegative 
real root). In this case, 

y7 2 + 4u(i? 0 - 1) b 
2 2 


r = 


( 5 ) 


where b = 


_ l-l _ fim) 

'Y+oi+fi Pfl+a/pfm+fi)) 


and u = 


(g+^V 

/3/3m(l+a/(7m+M)) 2 ' 
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Conditions involving the testing rate 

The rates 0 and a at which individuals are, respectively, deemed infected but still susceptible 
and diagnosed infected both depend on a treatment rate r. These rates will then depend on 
the testing rate and the test specificity and sensitivity. By applying the rates a = rip and 
6 = r( 1 — e) depending on the rate r, and the specificity e and sensitivity ip, we rewrite Rq: 

R = (P m r(l - e) + /3p)(rip + 7 m + p) 

° (r(l - e) + p)(rip + ^ + /x)( 7 m + p)' 

It is interesting to note that if r —> oo, Rq —> ( 7 ^+mu)4<{ i-e) = This expression 

is the equal to the Rq obtained for a SIR system at an infection rate (3 m and recovery rate 
j m . When r = 0, we have a regular SIR model formulated with parameters /?, 7 and p, 
such that Rq = 7 ^ 7 - Therefore, these are the two extreme values for Rq when varying the 
testing rate r. 

Conditions for R 0 < 1 

We consider now conditions for the disease not to go to the endemic state, he., for which 
Rq < 1. This condition would require that a polinomial p{r) < 0, where 

p(r) = (1 - e)ip(/3 m - 7 m - p)r 2 + 

(1 - e)( 7 m + p){(3 m - 7 - p) + pip(/3 - 7 ™ - p)r + 

/47m + ^)(/3-7-^)- 

The typical situation for analyzing p(r) < 0 is that an epidemic occurs if no treatment 
is realized, i.e. in this case Rq > 1, when r = 0. Therefore, we would expect the normal 
SIR threshold condition /3 > 7 + p,. Let us first assume that /3 m < 7 m + /x, which is the 
condition for the polynomial to be concave down. As shown when r —» 00 , this condition 
permits an Rq < 1 , which suggests that under treatment it is possible to bring the system 
to Rq = 1 for a finite r. Indeed, p(r) has two real roots, a positive root r c and a negative 
one. Therefore we would need r > r c in order to have p(r) < 0, and consequently Rq < 1. 

If we consider instead an asymptotic condition / 3 m /( 7 m + /x) > 1, given the reasonable 
assumptions > 7 and (3 > /3 m , would guarantee {3 > 7 + /x, and, by consequence, 

/ 3 >/ 3 m > 7 m + /x >7 + /x 

In this case, either there are no purely real roots or none of the roots are positive. Hence, 
p{r) > 0, Rq > 1. This essentially means that increasing the testing rate permits to decrease 
Rq, but not to avoid an epidemic. 

Different scenarios: numerical simulations 

Here we consider the model as given by Table [l] Values for parameters (3, /3 m , 7 , 7 m , /x 
were chosen such that we have Rq on a domain that permits situations from epidemic to 
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Figure 1. Countour plot under different testing rates. We show here testing rates 
from top to bottom: r = 0.02, r = 0.04, r = 0.06, r = 0.08. The intensity of the color 
indicates the Rq: colors closest to white are the smallest Rq values, whereas colors close to 
red indicate high values for Rq. 


controlled disease. We should note that under no treatment, we have an Rq condition: Rq = 

/3m 


= 1.456. When we have the limiting condition, r 


oo, we find Rq = = 0.654. 


Therefore these two values are the extreme points that we can find by considering different 
scenarios for 9 and a (or alternatively different values for r, e, and tjj). 

We first consider varying the sensitivity and specificity for different testing rates given 
by r. Figure [l] shows us from top to bottom the different countour plots for values of the 
basic reproduction number Rq when varying the testing rate r from r = 0.02 to r = 0.08. 
The intensity of the color permits us to observe the Rq values, in which a lighter-intensity 
color depicts smaller values for Rq. As expected if we increase the testing rate r, the plots 
become lighter, meaning that epidemics are less likely under higher testing rates. 

We also consider the cost ratio given by Eq. [2] and varying the testing rate r. Here, 
we consider the sensitivity equal to 0.98 and the specificity at 0.67, values found for a TB 


test by pooled results (Walusimbi et al. 2013), even though the model might be simple for 
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Figure 2. Cost ratio. Each of the curves depict the cost ratio obtained over time for 
different testing rates r. Curves vary from r = 0 (dark blue) to r = 10 (light blue). 

analyzing complex scenarios that TB presents. Figure [2] shows multiple curves, in which 
we observe for higher testing rates (lighter colors) the cost ratio increasing to higher peaks. 
This is expected because more individuals are under treatment. 

Finally, we show for the same testing rates as in Figure [2] how the epidemic develops 
(or not). We observe in Figure [3] that the incidence effectively decreases as the testing rate 
increases (lighter colors), which is the opposite effect of the cost ratio. Therefore, the ideal 
point should be a testing rate r that enables a basic reproduction number Rq < 1, but 
probably close to the threshold such that the cost ratio is not prohibitive. 

4 Discussion 

In this paper we establish a mathematical modeling that permits us to analyze the effects 
of imperfect testing on the dynamics of diseases. The modeling considers specificity and 
sensitivity of the testing. Therefore, given these two parameters we can find Rq, which tells 
us if the disease-free state is stable. In this case, even if few disease cases appear, we would 
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Figure 3. Incidence at different testing rates. Each of the curves show the number 
of infected persons under different testing rates varying from r = 0 (dark blue) to r = 10 
(light blue). The number of infected persons are found by summing variables I and I m . 


not expect an epidemic, to occur. Alternatively, we can find a testing rate that can bring 
the disease to a controled level for a specific test, for which specificity and sensitivity has 
been estimated. 

The testing rate, however, also raises the total treatment costs given a time interval. 
Hence, as expected, it is not viable to raise the testing rate indefinitely, for cost-effectiveness, 
especially when the per-capita treatment is already costly. 

We believe our model is particularly helpful to give insights about regions (controlled 
disease/epidemic, cost prohibitive/viable costs). We expect to extend the work to other 
similar models, in particular a SEIR model. This would be the case for diseases such as 
tuberculosis, and especially multi-resistant tuberculosis, which requires expensive treatment. 

We considered estimated rates for sensitivity and specificity typical of tests, such as a 
tuberculosis test GenExpert (Walusimbi et al. 2013), to evaluate possible scenarios. We 


also intend to study data from public health systems that show the impact of introducing 
treatment and the cost associated with the treatment. 
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The imperfect testing might bring individuals to states, for instance a false-positive 
infected, in which they might be diagnosed after a slow test and should return to susceptible 
state. We intend to treat this case in future work. 

We consider the expected cases in which /3 > j3 m , a smaller force of infection under 
treatment, and > 7 , a speedy recover rate under treatment. Those are reasonable 
assumptions in most cases. Different cases, however, can be considered for more complex 
scenarios. In the case of sexual-contact diseases, maybe some individuals might incur on 
reckless behavior, which would instead not increase the treatment, but increase the force 
of infection, as discussed by Wilson et al. (2008). We intend to pursue investigating such 
“anomalies” in the framework of our model. 
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